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1  AMTnACT  tUumim  too  — ^ 

ADC'i  Computer-Aided  Systems  Enginecriiig  Branch,  RBES,  has  .  jmented  on  in-housc 
effort  to  evaluate  the  structural  reliability  of  the  output  unveguido  window  on  a 
‘Traveling  Wave  Tube  (TWT).  This  window  acts  as  a  seal  between  the  TWT's  vacuum 
envelope  and  output  waveguide.  Its  purpose  is  to  prevent  any  loss  due  to  leakage  of 
the  vacuum  while  allowing  passage  of  the  microwave  signal.  This  particular  disk- 
shaped  window  is  constructed  of  a  ceramic  material,  bcrj^llj^  and  contains  an  inner 
ring  of  copper  and  an  outer  ring  of  Mon  1  K-500  (70-3pfNi-Cu) .  It  was  suspected  that 
excessive  thermal  stresses  associated  with  the  very  hi^  operating  tcraperatures  by 
this  window  had  caused  it  to  faiL< 

inite  element  analyses,  along  with  material  failure  theories  were  used  to  determine 
the  window's  response  to  a  time-dependent  heat  source  and  operating  heat  sink 
temperature^  Analyses  were  made  for  various  rates  of  heat  dissipation  through  the 
-window-ranging  from  20  to  120  watts.  It  was  determined  that  cracking  of  the 
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bcryllia  at  ics  outer  edge  could  occur  iC  more  than  20  watts  of  heat  were  to  be 
dissipated  through  the  window.  The  effect  of  power  cycling  on  the  window's 
ductile  materials  (copper  and  Monel)  was  also  studied.  It  was  found  that, 
depending  on  the  amount  of  heat  dissipation,  the  copper  portion  of  the  window 
could  survive  between  500  and  9000  power  cycles  before  failure  would  occur. 


PREFACE 


The  structural  reliability  of  a  ceramic  RF  output  window  of  a  Travel¬ 
ing  Wave  Tube  (TWT)  is  assessed  in  this  report.  This  TWT  was  part  of  the 
ALQ-99  ECM  system  currently  in  service  on  the  EF-in  aircraft.  The  Finite 
Element  program,  Numerically  Integrated  Elements  For  Systems  Analysis 
(NISA),  along  with  material  failure  theories  wore  utilised  for  this  anal¬ 
ysis.  All  of  the  work  was  accomplished  in-house  at  RADC. 

The  author  wishes  to  thank  Mr.  Edward  J.  Jones,  Ms.  Grotchen  A. 
Bivens,  Mr.  Douglas  J.  Holshauer,  and  Mr.  William  Bocchi,  RADC/RBES,  along 
with  Dr.  Edward  Danissewski,  RADC/OCTP,  for  their  technical  oasistonca  on 
this  project. 
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0.0  EXECUTIVE  SUHMAnX 


Tho  objective  or  thla  study  was  to  investigate  the  physical  integrity 
of  a  TWT's  output  RF  window  (Figure  1).  This  window  is  part  of  a  TVT  that 
is  used  in  tho  ALQ-99  ECM  system.  The  output  energy  that  this  window  is 
subjo'*''  to  causes  the  window  to  experience  very  high  operating  tempera¬ 
tures  that  cause  high  material  stresses  which  adversely  affect  its  reli¬ 
ability.  The  window  analyzed  in  this  study  had  been  taken  from  a  failed 
TVrr  and  was  found  to  be  cracked.  Many  mici'owave  tube  failures  have  been 
attributed  to  mechanical  stresses  induced  in  the  output  window  by  heat 
sources  which  induce  a  temperature  distribution  throughout  the  window. 
Any  quantitative  analysis  of  this  phenomena  Involves  the  calculation  of 
temperature  and  stress  distributions.  Heating  conditions  that  caused 
cracking  of  the  window  and  also  the  e'.*:'.'ct3  of  power  cycling  on  the  system 
(number  of  power  cycles  to  failure)  were  determined.  These  are  good 
indications  of  the  window's  design  reliability. 

In  order  to  assess  the  window's  reliability,  finite  element  simula¬ 
tions  were  employed.  Finite  Element  Analysis  (FEA)  is  a  computer  simula¬ 
tion  technique  which  predicts  the  physical  behavior  of  a  system  under  any 
specified  type(s)  of  load(3).  A  geometric  model  of  the  system  is  generated 
and  the  material  properties  of  the  model  are  defined  within  the  input  data 
file.  The  data  file  also  contains  thermal  loads  and/or  physical  con¬ 
straints  on  the  system  (depending  on  the  type  of  analysis  being  performed). 
The  data  file  is  then  read  by  the  finite  element  code  and  the  analysis  is 
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performed.  Several  parametric  studies  were  made  of  this  window  using  this 
technique  and  the  results  used  to  assess  the  window’s  reliability. 

The  results  of  the  parametric  studies  made  by  finite  element  analyses 
were  used  in  conjunction  with  the  Manson-Coffin  equation  to  predict  the 
number  of  power  cycles  the  copper  portion  of  the  window  would  survive  in 
the  ALQ-99  environment.  The  Manson-Coffin  equation  is  an  empirical  rela¬ 
tion  based  on  fatigue-life  data  which  considers  the  total  strain  within  a 
ductile  material  and  the  effect  of  load  reversals  on  these  strains.  One 
power  cycle  was  considered  to  represent  one  loading  cycle.  Therefore,  the 
number  of  power  cycles  to  failure  indicates  how  many  times  the  TWT  can  be 
operated  before  window  failure  would  occur. 

A  plot  was  obtained  relating  heat  dissipation  through  the  window  to 
power  cyc3  es  to  failure .  From  this  plot ,  for  any  value  of  heat  dissipation 
rate  (ranging  from  20  to  120  watts)  the  number  of  power  cycles  to  failure 
could  bo  determined  for  the  copper  portion  of  the  window.  Within  this 
heating  range,  the  number  of  power  cycles  to  failure  was  found  to  vary  from 
500  to  9000  cycles.  FEA  results  were  also  used  in  conjunction  with  the 
maximum  normal  stress  theory  for  brittle  fracture  to  determine  what  condi¬ 
tions  would  cause  the  beryllia  portion  cf  the  window  to  crack.  It  was 
determined  that  for  heat  dissipation  rates  through  the  window  of  greater 
than  20  watts,  the  beryllia  would  begin  to  crack  at  its  outer  circumfer¬ 
ence. 
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AH  cqnoluaiona  were  made  on  the  assumption  that  there  were  no  manu¬ 
facturing  defents  associated  with  the  window.  Failure  could  only  be  caused 
by  design  defects  or  tube  operation  beyond  design  capabilities.  Struc¬ 
tural  failure  is  but  one  failure  mode  associated  with  TWT's.  A  study  is 
currently  underway  at  RADC  to  identify  other  failure  modes  associated  with 
TWT's.  It  is  anticipated  that  the  results  of  these  studies  will  provide 
information  that  would  lead  to  design  and  manufacturing  of  more  reliable 
Traveling  Wave  Tubes. 
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1.0  INTRODUCTION 


A  Traveling  Wave  Tube  (TWT)  is  an  electron  tube  used  for  amplification 
of  microwave  and  millimeter  wave  signals.  Its  operation  depends  on  the 
interaction  of  a  beam  of  electrons  with  an  electromagnetic  wave.  Some 
applications  of  TWT's  include  use  in  space  communications,  electronic  war¬ 
fare,  radar  and  the  relaying  of  home  video  signals. 

Investigation  of  the  effect  of  high  operating  temperatures  on  the 
structural  rel^i ability  of  a  TWT's  output  window  was  the  main  objective  of 
this  study.  If  temperature  limits  are  exceeded  due  to  high  heat  dissipa¬ 
tion  and/or  inadequate  cooling  met}R.'dn.,  the  resulting  thermal  stresses  can 
reach  a  high  enough  value  to  cause  cracki tg  or  fatigue  failure  of  the 
window. 

The  most  widely  used  waveguide  window  for  high-power  tubes  is  the 
circular  ceramic  disk  window  in  a  section  of  round  waveguide.  This  type  of 
window  can  achieve  a  bandwidth  in  the  thirty  percent  range  and  has  a  high 
average  and  peak  power  handling  capability  (8).  This  ceramic  window  is 
constructed  of  beryllia  (BeO).  Bcryllia  ceramic  windows  are  used  for  high 
average  power  outputs  because  of  their  good  thermal  conductivity  at  high 
temperatures  (111  W/M-C  at  325°C  or  6^1.111  BTU/hr-ft-F  at  600°F) .  Beryllia 
is  avoided  in  many  cases  for  window  applications  because  of  its  toxicity, 
unless  its  high  thermal  conductivity  is  an  absolute  requirement.  The 
material  is  used  more  as  an  electrical  insulator  where  the  chance  of 
exposure  to  beryllia  in  powder  form  is  very  slim  (8). 
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The  system  analyzed  in  this  study  and  its  dimensions  arc  shown  in 
Figure  2.  This  system  consists  of  the  beryllia  disk  which  surrounds  a 
thin-walled  copper  tube  through  which  the  output  RF  energy  is  conducted. 
The  beryllia  disk  is  surrounded  by  a  nickel-copper  alloy.  A  laboratory 
analysis  determined  the  composition  of  this  alloy  to  be  nickel  and  36X 
copper.  This  material  was  assigned  the  properties  of  Monel  K-500  (70-30 
Mi-Cu),  the  material  whose  composition  most  closely  resembled  the  alloy's 
composition. 

When  the  TVfT  is  powered  up,  current  flow  through  the  copper  interacts 
with  the  copper  atoms  causing  displacement  and  heat  generation.  The  system 
experiences  a  uniform  temperature  distribution  which  induces  a  thermal 
stress  distribution  throughout  the  disk.  If  these  stresses  reach  a  high 
enough  value,  the  beryllia  disk,  which  behaves  in  a  brittle  manner,  can  be 
expected  to  crack. 

The  effect  of  thermal  cycling  from  repeatedly  powering  up  and  power¬ 
ing  down  the  system  can  lead  to  fatigue  failure  of  the  window.  When  the 
system  is  powered  up,  the  ductile  materials  in  the  window  (copper  and 
monel)  become  strained  due  to  thermal  stresses.  After  the  system  is 
powered  down,  the  materials  return  to  their  original  thermal  state.  How¬ 
ever,  not  all  of  the  total  strain  is  recovered.  A  portion  remains  within 
the  material  due  to  plastic  deformation.  Every  time  this  process  is 
repeated,  strain  accumulates  within  the  material  and  after  a  certain 
number  of  power  cycles,  the  total  strain  within  the  material  can  reach  a 
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great  enough  value  to  cause  failure.  The  predicted  number  of  power  cycles 
to  failure  gives  a  good  indication  of  the  window's  reliability. 

MIL-HDBK-ai7,  Reliability  Prediction  of  Military  Equipmen  ,  "ides 
a  model  to  predict  the  operating  failure  rate  for  Traveling  Uavi  Tubes 
based  on  tube  characteristics  and  operating  environment.  This  model  was 
formulated  from  failure  data  obtained  from  previously  designed  tubes  and 
only  considers  the  reliability  of  the  overall  THT.  Since  no  mathematical 
model  exists  that  can  be  used  to  predict  the  reliability  of  an  output 
window,  an  alternate  method  is  needed.  Therefore,  several  parametric 
studies  were  made  of  the  window  using  finite  element  techniques  and  the 
results  used  to  assess  the  window's  reliability.  Research  into  this  area 
turned  up  no  reports  of  previous  finite  element  analyses  done  on  TWT  output 
windows . 

The  FEA  code  used  at  RADC  is  Numerically  Integrated  Elements  for 
System  Analysis  (NISA).  In  this  study,  NISA's  transient  heat  transfer  code 
was  used  to  simulate  the  system's  response  to  a  heat  input  and  time  varying 
operating  heat  sink  temperature.  The  resulting  temperature  distribution 
was  used  as  a  thermal  load  into  NISA's  nonlinear  static  stress  module. 
Nonlinear  stress  analyses  were  performed  to  account  for  the  plastic  behav¬ 
ior  experienced  by  the  ductile  materials  of  the  system  beyond  their  respec¬ 
tive  yield  points.  Stresses  beyond  a  material's  yield  point  will  cause  the 
material  to  deform  in  a  plastic  manner.  Results  from  the  nonlinear  anal¬ 
yses  were  used  in  the  maximum  normal  stress  formula  to  predict 


8 


a  heating  value  that  would  cause  br^ittle  failure  of  the  beryllla  to  occur. 
These  results  were  also  used  in  the  Manson-Coffln  equation  in  order  to 
predict  a  number  of  power  cycles  to  failure  for  the  copper  portion  of  the 
window. 

This  report  is  organized  in  the  following  manner.  First  ,  a  de&‘  *p- 
tlon  of  the  finite  element  model  is  given.  Second,  the  inputs  iuto 
transient  heat  transfer  code  are  described  along  with  the  resultinj.  teia- 
perature  distributions  for  different  heating  values.  These  results  are 
then  input  to  the  nonlinear  code  which  gives  the  resulting  stress  distribu¬ 
tions  throughout  the  window.  Finally,  results  from  the  stress  analyses  are 
interpreted  using  material  failure  theories  for  both  brittle  and  ductile 
materials . 
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2.0  FINITE  ELEKEJfT  KODEL 


Finite  aXeacnt  wodcls  are  geometric  models  that  represent  an  actu?,l 
physical  system.  Since  heat  was  being  generated  on  the  entire  inner 
diameter  of  the  copper,  the  heat  transfer  was  axisymmotric;  that  is,  the 
temperature  varied  only  in  the  radial  direction  (x-direction)  (Figure  2). 
Because  the  window  is  of  uniform  thickness,  it  was  also  assumed  that  the 
temperature  remained  constant  throughout  the  thickness  of  the  disk  (y- 
direction).  With  these  assamptions,  the  heat  transfer  is  actually  one¬ 
dimensional,  however,  to  better  visualise  the  temperature  distribution 
throughout  the  system,  a  two-dimensional  model  was  created.  The  geometric 
model  is  illustrated  in  Figure  3.  Figure  3  is  ©."iented  such  that  the  x- 
axis  represents  the  radial  direction  and  the  y  axis  represents  the  axis  of 
rotation.  This  model  is  a  two-dimensional  cross-section  of  Figure  2  taken 
parallel  to  the  x-y  plane  (looking  into  the  z  axis). 

Using  HISA's  Geometry  Data  Base  Computer  Program,  this  model  was 
easily  created.  The  model  was  divided  into  four  patches  (regions)  with 
each  patch  representing  a  different  material.  Patch  1  represents  a  very 
elastic  material  used  to  simulate  boundary  conditions.  This  patch  was 
assigned  highly  elastic  properties  and  its  purpose  was  to  represent  the 
space  between  the  y-axis  and  the  copper  (patch  2) .  The  space  was  repre¬ 
sented  in  this  way  so  that  the  model  could  be  properly  constrained  in  the  x 
and  y  directions.  Constraining  the  left-hand  side  of  the  copper  would  have 
been  invalid  since  the  whole  model  must  be  allowed  to  expand  radially. 


10 


K 


11 


FIGURE 


tv.'xpPi'Uv^ly. 

The  finite  elraent  nedel  (FiRure  H)  in  defin<'d  by  twa  isabheK-itleal 
p3ra:3eters,  n^'ins  and  elenenta.  A  nede  la  a  Rrid  p^lnt  l^seated  in  npaee. 
Pinco  this  KOvlol  was  created  in  t«d  diaensicos,  each  nc-de  was  defined  by 
twa  caardinates,  x  and  y.  The  eleaent  is  defined  by  a  material  index 
nuabe  and  a  sot  af  nodes.  Each  index  number  has  specific  material  proper¬ 
ties  ssoeiated  with  it.  The  p:*opertie3  can  be  defined  to  be  constant  or 
temperature  dependent.  Two  of  the  materials,  e.  pper  and  boryMia,  had 
temperature-dependent  thermal  eonductivities.  A  few  sets  of  data  polnt.s 
were  input  into  the  data  file  and  the  NISA  code  linearly  interpolated 
between  these  points  to  determine  the  conductivity  at  any  temperature. 
Only  one  value  of  thermal  conductivity  could  bo  found  for  monel,  therefore, 
its  thermal  conductivity  was  assumed  to  be  independent  of  temperature. 
Since  the  thermal  conductivities  of  most  metals  do  not  depend  heavily  on 
temperature,  this  is  a  reasonable  assamption. 

The  elements  used  in  this  study  wci*o  two-dimensional,  four-sided, 
axisynmetric  solid  elements.  Each  element  was  defined  by  four  nodes,  one 
at  each  corner.  Two  types  of  these  elements  were  utilized,  the  first  for 
the  heat  transfer  analysis  and  the  second  for  the  nonlinear  static  stress 
analysis.  These  elements  were  identical  in  shape  and  orientation  with 
their  only  difference  being  that  the  heat  transfer  element  has  only  one 
degree  of  freedom,  temperature,  while  the  static  element  has  two  degrees  of 
freedom,  movement  in  both  the  x  and  y  directions. 
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Finite  Element  Model 
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R6)«.^ns  That  Represent  The  Materials  Of  The  Model 


FIGURE 
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The  finite  eie:aent  aodcl  containe  a  tatal  of  800  nodes  and  690  cle- 
Rents.  Bach  aateriaX  contains  the  following  node  and  olcaent  nuxbers: 


Material 

Node 

Element  i's 

Copper 

1-6»l 

1-95 

Beryllia 

65-980 

96-920 

.Monel 

981-736 

921-695 

Elastic 

737-800 

696-690 

The  results  obtained  for  the  clastic  aaterial  wore  reviewed  to  verify 
the  accuracy  of  this  technique,  but  not  considered  in  the  stress  analyses. 

For  both  heat  transfer  and  static  stress  analyses,  all  nodes  lying  on 
a  boundary  between  materials  were  coupled  to  each  other.  What  this  means 
is  that  for  the  thermal  analyses,  those  nodes  were  required  to  have  thn 
same  temperature  and  for  the  stress  analyses,  they  wore  requi.'^d  ta  move 
equal  amounts  in  both  the  x  and  y  directions.  If  these  nodes  were  merged 
instead  of  coupled,  the  temperature  and  stress  results  given  would  have 
been  average  values  between  nodes  of  two  dissimilar  materials,  which  would 
have  been  inaccurate.  If  these  nodes  were  left  uncoupled,  the  HISA  pro¬ 
grams  would  not  execute  due  to  disjointed  elements  at  the  bounda.''io3. 


3.0  TRANSIE^IT  HEAT  TRANSFER  ANALTSIS 


Model 


A  transient  heat  transfer  analysis  simulates  the  model's  thermal 
statvT  when  the  temperature  distribution  varins  with  time.  Since  the  window 
is  enclosed  in  a  section  of  round  waveguide,  it  was  assumed  that  the 
primary  mode  of  heat  transfer  was  conduction.  There  was  no  convection  or 
radiation  to  the  outside  air.  The  thermal  resistance  associated  with 
conduction  is  considerably  less  than  the  resistances  for  convection  and 
radiation,  therefore,  this  is  a  reasonable  assumption.  This  case  of  pure 
conduction  represents  a  worse  case  scenario  under  which  the  window’s  reli¬ 
ability  will  be  characterized. 

This  analysis  was  performed  with  heat  being  generated  on  the  entire 
inner  diameter  of  the  copper.  The  heat  generation  was  simulated  as  concen¬ 
trated  nodal  heat  fluxes  at  nodes  1,  5,  9.  13,  17,  21,  29,  33,  37,  Ml,  *15, 
53  and  6l  (see  Figure  *1).  In  order  to  provide  a  path  for  the  heat  to  flow, 
a  heat  sink  temperature  was  designated  at  the  outer  edge  of  the  monel. 
These  nodes  represent  the  extreme  outer  edge  of  the  model  (see  Figure  *J). 
The  heat  sink  represents  an  area  where  the  system  is  cooled  to  a  specified 
temperature  to  prevent  possible  melting  of  any  components.  A  tube  operat¬ 
ing  temperature  profile  was  measured  during  an  actual  ALQ-99  mission 
flight  test  of  the  tube  and  is  shown  in  the  curve  of  Figure  5.  This  data 
was  obtained  from  temperature  readings  taken  on  the  top  of  TWT's  beam 
collector.  The  location  of  the  collector  is  close  to  the  outer  edge  of  the 
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Environmental  Flight  Test  Profile 


FIGURE 


SF  window,  therefore,  it  was  assumed  that  the  system's  heat  sink  tempera¬ 
ture  followed  this  time-amplitude  curve.  A  number  of  data  points  were 
inserted  into  the  input  file  to  represent  this  curve.  The  curve  was  then 
referenced  by  the  sink  temperature  data  group  in  the  NISA  input  f  le  to 
give  temperature  and  heat  source  values  at  any  time. 

Results 


NISA's  transient  heat  transfer  program  was  executed  in  time  steps  of 
ten  seconds  over  a  span  of  580  seconds.  Temperature  contours  for  50  and 
100  watt  heat  sources  are  shown  in  Figures  6a  and  6b  (all  temperature  units 
are  given  in  degrees  C).  These  plots  show  the  system's  thermal  state  at 
the  conclusion  of  the  run  (580  seconds).  It  is  apparent  that  the  widths  of 
the  temperature  bands  through  the  beryllia  are  wider  than  through  the  other 
materials.  This  is  because  of  beryllia' s  much  higher  thermal  conductiv¬ 
ity.  Figures  7a  and  7b,  obtained  through  NISA's  post  processor,  are  time 
vs.  temperature  curves  at  four  nodes  in  the  model  for  heat  dissipation 
rates  of  50  and  100  watts  respectively.  Temperature  values  at  nodes  1,  65 
and  i|8l  were  plotted  because  these  are  the  nodes  of  maximum  temperature  for 
the  copper,  beryllia  and  monel,  respectively.  In  all  cases,  the  tempera¬ 
ture  of  the  model  reaches  a  maximum  at  approximately  290  seconds  and 
remains  relatively  constant  for  the  duration  of  the  run.  Node  1  is  the 
node  of  greatest  temperature  change  from  the  heat  sink  since  it  is  one  of 
the  nodes  where  heat  is  being  generated.  Nodal  temperatures  decrease  with 
increasing  distance  from  the  heat  source. 
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Th«  results  of  these  heat  transfer  analyses  have  shown  that  at  any 
titts  during  the  TWT's  operation  i  the  teaperature  distribution  throu?,hoot 
the  Sr  window  is  in  the  radial  direction  and  does  not  vary  in  the  vertical 
direction. 
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Temperature  Distribution  -  50  Watts 


G3  G3 


U" 

a 

g 

o 

c* 


nj  ni 

■  CD  Q 

*  * 

H  UJ  UJ 
a  US  V 
lu  nji  r- 
X  •  • 

♦-  (U 


cr 

Qc 


iimiiiiiiiiii 

iimiiiiieiiii 

IIIMIIIIIIIIOI 


ifiiiiiiiiiiiii 


CO 

(0 


III 


O 


11. 


19 


125.6 

Erf?C  -HISfi^DISPLAY 


ISOTHCRM  CONTOURS 

Temperature  Distribution-100  Watts  transient  heat 
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Time  vs  Temperature  Curves  50  Watts 
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FIGURE 


Time  vs  Temperature  Curves  100  Watts 


»"  M  C3  04 
Q  ?a  SO  o 

UJ  »  •  »  * 

lU  UJ  U  UJ 
CS  O  to  so  G3 

a  lA  fU  ®  03 

yi  fu  T  CD  lA 

M 

X 


u  Q,  O  ui  O 

Z  C  r  E  r 

»-i  uj  a:  M  «r 

h-  F  os  H  a: 


-  lA 
so 


^  hi 

tn 


(D 

S) 

(S 

o 

(S 

• 

• 

• 

• 

S) 

1 

in 

in 

in 

in 

in 

ru 

't 

so 

00 

GO 

in 

'T 

cn 

OJ 

T— 

r— 

OJ 

o  930)  3ani«a3cjW3i  lyoow 


22 


FIGURE 


»<.0  THERMAL  STRESS  ANALYSES 


Thermal  stresses  arc  stresses  arising  from  temperature  variations 
within  a  material.  When  a  material  is  subjected  to  a  temperature  gradient, 
the  various  fibers  tend  to  expand  different  amounts.  This  expansion 
induces  tensile  and  compressive  stresses  throughout,  if  the  material  is 
physically  constrained.  When  these  stresses  exceed  the  yield  strength, 
deformation  and/or  failure  will  occur. 

The  results  of  NISA's  transient  heat  transfer  analysis  gave  tempera¬ 
tures  at  each  node  for  a  given  heat  source.  The  resulting  output  file  was 
saved  at  the  time  step  where  maximum  temperature  occurs  for  each  source  and 
used  as  a  thermal  loading  condition  into  NISA's  nonlinear  static  stress 
module.  This  time  step  was  chosen  at  the  conclusion  of  each  run  (580 
seconds)  which  also  corresponded  to  the  steady  state  condition.  These 
nonlinear  stress  analyses  take  into  account  the  behavior  of  a  ductile 
material's  stress-strain  curve  beyond  its  yield  point.  Data  points  con¬ 
taining  stress  and  corresponding  strain  values  for  both  copper  and  monel 
beyond  their  respective  yield  points  were  inserted  into  the  input  file. 
Copper's  stress-strain  curve  beyond  its  yield  point  is  shown  in  Figure  12. 
The  same  time-amplitude  curve  used  to  represent  the  sink  temperature  in  the 
transient  heat  transfer  analyses  was  used  to  represent  the  temperature 
loading  in  these  nonlinear  stress  analyses. 

A  series  of  thermal  stress  analyses  were  performed  on  the  RF  output 
window.  First,  stress  analyses  of  the  window  were  made  when  the  entire 
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oystea  was  assigned  a  finite  tcaporaturo  (zero  heat  generation).  This  was 
done  to  dotcralne  if  craoking  of  the  beryllia  could  possibly  occur  when  the 
TWT  was  in  the  power  off  condition  and  subject  to  extremely  cold  ambient 
temperatures.  K  nonlinear  static  stress  analysis;  was  made  with  the  entire 
temperature  of  the  system  being  dropped  to  -150®C  (approximately  -83®F). 
The  results  showed  the  maximum  principal  stress  to  be  8i4.86N/sq  mm  and  the 
minimum  principal  stress  to  be  -192.33N/3q  ram.  Since  neither  of  these 
values  approach  the  tensile  or  compressive  stress  of  beryllia  (159  and 
-130,000N/3q  mm,  respectively),  the  beryllia  will  not  crack  under  this 
condition  X 

Next,  stress  analyses  were  made  using  the  heat  transfer  results  for 
heat  dissipation  rates  ranging  from  20  bo  120  watts.  These  results  allowed 
several  plots  to  be  developed  relating  thermal  stresses  to  elapsed  time 
after  power  up  of  the  system.  Because  of  a  lack  of  information  regarding 
stress-free  temperature  for  these  materials,  it  was  assumed  that  a  state  of 
zero  stress  existed  at  25°C  (approximately  room  temperature).  Stresses  at 
each  node  were  dependent  upon  the  difference  between  the  temperature  at 
that  node  and  this  stress-free  temperature.  The  results  do  not  consider 
the  residual  stresses  that  exist  prior  to  power  up.  Residual  stresses,  in 
general,  are  caused  by  prior  thermal  cycling,  manufacturing  processes, 
creep,  etc.  Since  copper  and  monel  were  assumed  to  behave  as  ductile 
materials,  and  beryllia  a  brittle  material,  different  methods  of  determin¬ 
ing  the  possibility  of  failure  had  to  be  employed. 
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Bcpyllla 


Failure  of  bribble  maberials  occura  due  bo  a  peak  load.  These  mabe- 
rials  do  nob  fall  due  bo  accumulablon  of  fablgue,  as  do  duoblle  mabtirials. 
To  predlcb  bribble  fraobure,  bhe  maximum  normal  sbress  bheory  was  used. 
This  sblpulabes  bhab  fraobure  in  a  bribble  maberial  will  occur  whenever  bhe 
maximum  principal  sbress  ab  a  poinb,  if  lb  is  in  benslon,  is  equal  bo  or 
greaber  bhan  bhe  bensile  sbrengbh  of  bhe  maberial.  Similarly,  if  bhe 
specimen  is  in  compression,  failure  will  occur  if  bhe  algebraic  minimum 
principal  sbress  ab  a  polnb  is  equal  bo  or  more  negablve  bhan  bk^  mabe- 
rial's  compressive  sbrengbh  (7).  Principal  sbresses  are  normal  sbresses 
assoclabed  wlbh  axes  on  whose  normal  planes  bhere  exlsb  exbreme  normal 
sbress  and  zero  shear  sbress.  Ib  is  generally  bhe  case,  and  bhis  case  is 
no  excepbion,  bhab  maberials  have  greaber  sbrengbh  in  compression  bhan  in 
benslon.  Beryllla  has  a  bensile  sbrengbh  of  159  N/sq  mm  and  a  compressive 
sbrengbh  of  130,000  N/sq  mm.  (Appendix  A  conbains  properbies  for  all 
maberials  used  in  bhis  sbudy.) 

Figures  Pa  and  8b  are  plobs  of  maximum  principal  sbress  (N/sq  mm) 
bhrough  bhe  beryllla  wibh  50  and  100  wabbs  of  heab  being  dissipabed  bhrough 
bhe  window,  respecbively.  All  posibive  sbress  values  indicabe  benslon, 
while  all  negabive  sbress  values  indicabe  compression.  These  con hours 
demonsbrabe  bhe  increased  sbress  levels  associabed  wibh  increased  heab 
dissipabion.  Higher  heab  dissipabion  rabes  induce  more  posibive  bensile 
sbresses  and  more  negabive  compressive  sbresses.  Figures  9a  and  9b  are 
curves  of  bime  vs.  maximum  principal  sbresses  ab  node  142  in  bhe  beryllia 
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ab  50  and  100  watts  heat  dissipation.  Hode  is  located  at  the  beryllla- 
monel  interface  and  is  the  node  of  maximum  principal  stress  (maximum  ten¬ 
sion)  for  the  beryllia.  Tensile  stresses  are  greatest  at  the  nodes  near 
the  outer  diameter  of  the  beryllia  because  as  the  system  is  heated,  the 
monel  undergoes  greater  expansion  than  does  the  beryllia  (this  is  because 
monel  has  a  greater  coefficient  of  thermal  expansion  than  beryllia),  caus¬ 
ing  the  beryllia  to  bo  in  tension.  When  heat  is  being  dissipated  through 
the  window  at  a  rate  of  20  watts,  a  maximum  principal  stress  of  156.5N/sq 

mm  was  calculated.  Since  this  value  approaches  beryllia' s  tensile 
2 

strength  of  159N/mm  ,  oracking  of  the  beryllia  can  be  expected  to  occur  for 
heat  dissipation  rates  above  20  watts. 
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STRESS  CQHTOURS 

Stress  Distribution  -  BeO  50  Watts  si  priucpl  stress 
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STRESS  CONTOURS 
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In  ord‘?r  U  |}j*«dlot  fchfl  lift  of  the  copper  partioo  of  tJie  BF  window, 
the  effect  of  fatigue  due  to  power  cycling  was  studied.  The  systeis  experi¬ 
ences  one  power  cycle  when  it  goes  f.'^M  its  initial,  or  stress-free  ten- 
perature,  to  its  power  up  tewperature,  then  back  to  its  original  state. 
If,  when  the  systen  is  powered  up,  the  stresses  thi'ough  the  copper  exceed 
its  yield  stress,  the  copper  will  experience  a  certain  amount  of  strain. 
After  the  aysteoi  is  returned  to  its  original  thensal  state,  only  the 
elastic  strain  is  recovered.  The  strain  that  regains  within  the  Material, 
the  plastic  strain,  aeeuMulates  with  every  power  cycle.  Therefore,  after  a 
certain  nunber  of  power  cycles,  the  strain  will  reach  a  high  enough  value 
to  cause  the  copper  to  fail. 

When  determining  the  total  strain,  the  following  procedure  is  used. 
First,  Von-Mises  stresses  in  the  copper  are  determined.  The  Von-Mises 
stress  is  a  combined  stress  value  that  is  widely  accepted  in  predicting 
failure  of  ductile  materials.  Von-Mises  stresses  through  the  copper  were 
calculated  for  different  vtlues  of  heat  dissipation.  Figures  10a  and  10b 
are  contour  plots  of  Von-Mlses  stresses  in  the  copper  for  50  and  100  watts 
of  heat  dissipation.  As  was  the  case  with  beryllia,  the  stress  levels 
increased  with  increased  heat  dissipation  through  the  window.  The  Von- 
Kises  stress  is  a  maxiaimi  at  nodes  <1  and  6b,  the  lower  and  upper  comers  of 
the  copper- beryllia  interface,  respectively.  Since  the  distribution  is 
syasaetrlc  about  the  center  of  the  model,  these  stress  values  are  Identical 
and  only  the  data  for  node  b  was  taken  into  account.  The  Von-Mlses 
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Stress  Distribution-Cu  100  Watts 
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FIGURE  10b 


stresses  are  plotted  in  Figures  11a  and  11b  against  elapsed  time  after 
power  up  at  node  <1.  This  curve  demonstrates  the  increase  in  the  Von-Mises 
stress  up  until  approximately  290  seconds,  when  the  TVTT  operating  tempera¬ 
ture  reaches  a  maximum  (see  TVTT  opei'ating  profile  in  Figure  5)  •  After  this 
time,  the  stresses  increase  only  sli^tly  as  time  is  increased. 

The  total  strain  is  then  obtained  for  a  given  stress  from  the  stress- 
strain  curve  for  copper  shown  in  Figure  12.  This  total  strain  value  is 
translated  into  number  of  cycles  to  failure  using  the  Manson-Coffin  rela¬ 
tion  for  a  specific  material.  The  Hanson  Coffin  model,  along  with  the 
parameter  values  for  copper,  is  shovfn  in  Figure  13.  These  values  were 
inserted  into  the  Manson-Coffin  model  and  a  plot,  shown  in  Figure  1*1  was 
generated.  This  plot  contains  two  lines  and  one  curve.  The  lines  repre¬ 
sent  the  elastic  and  plastic  components  and  the  curve  represents  the  sum  of 
these  two  components.  To  obtain  the  number  of  cycles  to  failure  for  a 
given  strain,  the  total  strain  value  at  node  *1  is  plotted  onto  the  curve 
and  the  nmaber  of  cycles  to  failure  is  determined.  The  chart  below  lists 
maximum  Von-Mises  stresses  (node  *1),  strains,  and  power  cycles  to  failure 


for  different  values 

of  heat  dissipation  rates. 

HEAT  (WATTS) 

STRESS  (N/sq  mm) 

STRAIN 

CYCLES 

20 

J19.18 

.0054 

9000 

50 

56.64 

.007 

1450 

80 

59.62 

.011 

950 

90 

61.49 

.0118 

800 

100 

63.17 

.0125 

700 

120 

65.96 

.0145 

500 
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Stress-Strain  Curve  For  Copper  (Above  Yield  Point) 
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Manson-Coffin  Model 


ELASTIC 

STRAIN 

PARAMETER 

‘f 

E 

b 


PLASTIC 

STRAIN 

COPPER 

345  N/sq  mm 
.3 

1.25  E6  N/«q  mm 
-.05 
-.6 


;  FATIGUE  STRENGTH  COEFFICIENT 
‘f  :  FATIGUE  DUaiLITY  COEFFICIENT 
E  :  MODULUS  OF  ELASTICITY 
b  :  FATIGUE  STRENGTH  EXPONENT 
c  :  FATIGUE  DUCTILITY  EXPONENT 
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Manson-Coffin  Curve  For  Copper 
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A  plot  of  heat  dlsaipatlofi  through  tho  window  (normal  scale)  vs. 
number  of  power  cycles  to  failure  (log  scale)  for  the  copper  is  shown  in 
Figure  15.  This  plot  indicates  that  an  increased  heat  input  will  cause  a 
decrease  in  the  number  of  cycles  to  failure,  indicating  a  loss  reliable 
system.  By  using  this  graph  for  a  given  heat  input,  a  namber  of  power 
cycles  to  failure  can  be  determined. 
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Cycles  To  Failure  vs  Heat  Dissipation  -Copper 
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Monftl 


Like  copper,  mcoel  behaves  in  a  ductile  manner.  Vwi-Miaea  atressca 
were  calculated  and  their  contours  are  plotted  in  Figures  l6a  and  l6b  for 
heat  sources  of  50  and  100  watts,  reapectively.  These  contour  shapes 
remain  approximately  the  same  for  all  cases  with  the  larger  stress  magni¬ 
tudes  occurring  for  larger  heat  sources. 

Time  vs.  stress  curves  wore  then  plotted  for  50  and  100  watt  heat 
sources  at  node  *I81  (Figures  17a  and  I7b).  The  largest  Von-Miscs  stresses 
through  the  monel  occurred  at  this  node  which  is  located  at  the  lower 
corner  of  the  beryllia-moncl  interface. 

Monel's  iield  stress  of  758N/sq  mm  was  not  exceeded  for  any  of  the 
nonlinear  statiu  runs  that  were  performed  (up  to  120  watts).  Therefore, 
this  portion  of  the  window  is  not  expected  to  contribute  to  the  RF  window's 
failure. 
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STRESS  CONTOURS 

Stress  Distribution  -  Monel  50  Watts  von-mises  stress 

VIEW:  9.19E  +  00 


Stress  Distribution  -  Monel  100  Watts  von-mises 

VIEW:  2.13E+01 

RANGE:  3.80E+02 
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Time  vs  Stress  Curve  -  Node  481  100  Watts 
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Mierowavn  tubfl  reliability  is  strongly  dependent  on  three  factors. 
First,  defects  Introduced  during  the  manufacturing  process  adversely 
affect  reliability.  Produolbillty  problems,  poor  workmanship  and  lack  of 
process  control  are  major  contributors  to  manufacturing  defects. 
Secondly,  tube  reliability  is  heavily  dependent  upon  operating  procedures 
and  handling.  Finally,  adequate  design  margin  must  exist  between  the 
operating  point  and  the  ultimate  design  capability  of  the  tube  in  order  to 
have  reliable  operation. 

At  the  onset  of  this  study,  it  was  assumed  that  neither  manufacturing 
processes  nor  handling  procedure  played  a  part  in  the  failure  of  this 
window.  It  was  ‘.he  purpose  of  this  study  to  investigate  the  Impact  of 
increased  operating  temperature  on  its  stru^ural  integrity. 

The  two  most  significant  conclusions  made  based  on  the  results  of 
these  analyses  were: 

-  At  heat  dissipation  rates  greater  than  20  watts,  cracking  of  the 
beryllia  would  begin  to  occur. 

•>  Pep*;*?  -  ■  mount  of  heat  being  dissipated  through  the  sys¬ 
tem,  the  ;  O'  *.  on  of  tne  window  could  be  expected  to  survive 

anywhere  between  IJOO  and  9000  power  cycles. 
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Some  other  conclusions  that  can  be  made  fram  these  results  arc: 


-  Cracking  of  the  beryllia  would  occur  in  tension,  there  was  no 
evidence  found  that  would  predict  compressive  failure. 

-  A  stress  analysis  of  the  system  in  the  power  off  condition  subject 

^  to  an  extremely  cold  ambient  temperature  (-150°C)  indicated  that  the 

beryllia  would  nob  fail  under  this  extreme. 

-  Stresses  through  the  monel  portion  of  the  window  did  not  approach 
its  yield  stress  for  any  of  the  analyses  made.  Therefore,  this  region 
will  not  contribute  to  the  window's  failure. 

The  results  of  this  study  have  demonstrated  that  finite  element  anal¬ 
ysis  is  a  valuable  tool  that  can  be  used  to  predict  structural  reliability 
of  this  RF  window.  Lack  of  sufficient  TWT  data  did  not  permit  a  conclusion 
to  be  made  on  whether  inadequate  design  or  operation  beyond  design  limits 
had  caused  the  window  to  fail.  It  should  be  noted  that  this  study  was 
/  limited  to  an  output  window  of  this  size,  composition  and  operating  envi¬ 

ronment  . 

Finite  element  analysis  can  prove  to  be  a  valuable  tool  in  advancing 
TWT  technolc 'V  by  aiding  in  the  design  of  future  TWT  components.  Up-front 
structural  reliability  assessments  can  be  made  on  a  component  before  pro¬ 
duction  actually  begins  to  insure  that  it  can  withstand  the  conditions  it 
was  designed  for.  Finite  element  analysis  can  also  be  used  as  an  analysis 
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APPENDIX  A 


MATERIAL  PROPERTIES 


COPPER 

DBR7LLIA 

Material  ID=1 

Material  ID=2 

B=l.^n  E+05 

E=3.25  E+05 

»U=.3R 

NUs  .23 

ALPHAS 17. 6  E-06 

ALPHAsS.O  E-06 

k=.389  e  -73  DEO  C 

TENS.  ST. =159.0 

.1101  e  27  "  " 

COMP.  ST.=1.3E+05 

.R09  e  127  "  " 

k=.272  e  27  DEG  C 

.«124  P  327  '•  ” 

.196  e  127  "  " 

7=  ui.n 

.111  §  327  " 

0=  365 

.07  e  527  "  " 

DEHSs  8.»17  E-06 

DEllSs  3.0E-06 

MONEL  K-500 

ELASTIC 

Material  10=3 

Material  IDs'l 

E^l.79  E+05 

E=2.76  E+0'1 

IIU=.32 

HUs.il 

ALPHAS 13. 68  E-06 

ALPHA=5.0  E-0'1 

7=758.5 

k=.017JI5 

C=  '118.7 

DEN3=  S.‘17  E-06 

S7MB0L  AND  DEFINITION 

UNITS 

E; 

Modulus  of  Elasticity 

N/sq  mm 

NU: 

Poissons  Ratio 

Dimensionless 

ALPHA: 

Coefficient  of  Thermal  Expansion 

1/DEG  C 

7; 

7ield  Strength 

N/sq  mm 

TENS  ST: 

Tensile  Strength 

N/sq  mm 

COMP  ST: 

Compressive  Strength 

N/sq  mm 

k: 

Thermal  Conductivity 

Watts/mm  C 

C: 

Specific  Heat 

Joules/kg  C 

DENS: 

Mass  Density 

kg/cubic  mm 
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RADC  plans  and  executes  research,  development,  test  and 
selected  acquisition  programs  in  support  of  Command,  Control, 
Communications  and  Intelligence  (C*I)  activities.  Technical  and 
engineering  support  un'thin  areas  of  competence  is  provided  to 
BSD  Program  Offices  (POs)  and  other  ESD  elements  to 
perform  effective  acquisition  of  C*/  systems.  The  areas  of 
technical  competence  include  communications,  command  and 
control,  battle  management  information  processing,  surveillance 
sensors,  intelligence  data  collection  and  handling,  solid  state 
sciences,  electromagnetics,  and  propagation,  and  electronic 
reliability /maintainability  and  compatibility. 


